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We summarize the top-quark mass measurements from the CDF and D0 experiments at 
Fermilab. We combine published Run I (1992-1996) measurements with the most recent 
preliminary Run II (2001-present) measurements using up to 2.8 fb _1 of data. Taking 
correlated uncertainties properly into account the resulting preliminary world average 
mass of the top quark is M t = 172.4±0.7 (stat.)ztl.O (syst.) GeV/c 2 , assuming Gaussian 
systematic uncertainties. Adding in quadrature yields a total uncertainty of 1.2 GeV/c 2 , 
corresponding to a relative precision of 0.7% on the top-quark mass. 



1 The Tevatron Electroweak Working Group can be contacted at tev-ewwg@fnal.gov. 
More information can be found at [EttpT77tevewwg.fnal.gov 
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1 Introduction 



The experiments CDF and D0, taking data at the Tevatron proton-antiproton collider lo- 
cated at the Fermi National Accelerator Laboratory, have made several direct experimental 
measurements of the top-quark pole mass, M t . The pioneering measurements were based on 
about 100 pb" 1 of Run I (1992-1996) data [Il[2l[3l[llEl[6l[Tl[8l[9l UHl [HI Ll2l and include 
results from the it — *> qq'bqq'b (all-j), the it —> ivqq'bb (1+j), and the it — > i+vbt'Vb (di- 
1) decay channels^]. The Run II measurements summarized here are the most recent results 
in the 1+j, di-1, and all-j channels using 1.9 — 2.8 fb™ 1 of data and improved analysis tech- 
niques [i3i [iii rrsi iisi nn na nn [20] . 



This note reports the world average top-quark mass obtained by combining five published 
Run I measurements [21 El El UHl EEE] with four preliminary Run II CDF results [131 CHI 
[T5l [T6] and three preliminary Run II D0 results [T71 HH1 HI [20]. The combination takes into 
account the statistical and systematic uncertainties and their correlations using the method of 
references [2D [22] and supersedes previous combinations [231 EH ESI ESI EH EE] • 

The input measurements and error categories used in the combination are detailed in Sec- 
tions E] and El respectively. The correlations used in the combination are discussed in Section @] 
and the resulting world average top-quark mass is given in Section EJ A summary and outlook 
are presented in Section El 



2 Input Measurements 



For this combination twelve measurements of M t are used: five published Run I results, and 
seven preliminary Run II results, all reported in Table [H In general, the Run I measurements 
all have relatively large statistical uncertainties and their systematic uncertainty is dominated 
by the total jet energy scale (JES) uncertainty. In Run II both CDF and D0 take advantage 
of the larger it samples available and employ new analysis techniques to reduce both these 
uncertainties. In particular the JES is constrained using an in-situ calibration based on the 
invariant mass of W — > qq' decays in the 1+j and all-j channels. The Run II D0 analysis in 
the 1+j channel constrains the response of light-quark jets using the in-situ W — > qq' decays. 
Residual JES uncertainties associated with 77 and pr dependencies as well as uncertainties 
specific to the response of 6-jets are treated separately. Similarly, the Run II CDF analyses 
in the 1+j and all-j channels also constrain the JES using the in-situ W — > qq' decays. Small 
residual JES uncertainties arising from rj and dependencies and the modeling of 6-jets are 
included in separate error categories. The Run II CDF and D0 di-1 measurements use a JES 

2 Here £ = e or /1. Decay channels with explicit tau lepton identification are presently under study and are 
not yet used for measurements of the top-quark mass. 



2 





Run I published 


Run II preliminary 






CDF 




D0 






CDF 






D0 






all-j 


1+j 


di-1 


1+j 


di-l 


1+j 


di-1 


all-j 


trk 


1+j /a 


1+j/b 


di-1 


JCdt 


0.1 


0.1 


0.1 


0.1 


0.1 


2.7 


1.9 


2.1 


1.9 


1.0 


1.2 


2.8 


Result 


186.0 


176.1 


167.4 


180.1 168.4 


172.2 


171.2 


176.9 


175.3 


171.5 


173.0 


174.4 


IJ-CjO 


(J. U 


U.U 


U.U 


0.0 


0.0 


U.y 


n n 
U.U 


1 
i.y 


U.U 


n 7 
U. / 


1 A 

1.4 


U.U 


o TT?Q 
aJUjO 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


n n 
U.U 


n n 
U.U 


0.0 


0.8 


0.8 




DJUjD 


O.D 


0.D 


0.8 


0.7 


0.7 


n A 

0.4 


U.4 


U.6 


0.0 


0.0 


0.1 


0.3 


CJrvo 


n 

0.0 


O 7 

2.7 


2.b 


2.0 


2.0 


n 

0.3 


1 7 
1. 1 


U.u 


U.6 


U.U 


0.0 


U.U 


UJUjO 


U.o 


U. 1 


U.O 


0.0 


0.0 


n n 
U.U 


u. 1 


n 1 

U.l 


U.U 


U.o 


n n 
U.U 


l.b 


t- TTTQ 


4.U 


0.4 


Z. 1 


2.5 


1.1 


n A 
U.4 


1 9 


6 


n 1 
U.l 


U.U 


n n 
U.U 


U.U 


1ot-iP+ 
iep-r 1 


n n 

u.u 


n n 

u.u 


n n 

u.u 


0.0 


0.0 


n 9 


0.1 


0.0 


1 1 

± . X 


n n 

U.U 


n n 

U.U 


n n 

U.U 


>I 1 p*ri a ] 


1.8 


2.6 


2.8 


1.1 


1.8 


0.3 


0.8 


0.5 


1.6 


0.5 


0.5 


0.5 


BG 


1.7 


1.3 


0.3 


1.0 


1.1 


0.4 


0.4 


1.0 


1.6 


0.4 


0.4 


0.6 


Fit 


0.6 


0.0 


0.7 


0.6 


1.1 


0.2 


0.6 


0.5 


1.4 


0.3 


0.2 


0.3 


MC 


0.8 


0.1 


0.6 


0.0 


0.0 


0.5 


0.9 


0.5 


0.6 


0.0 


0.0 


0.0 


UN/MI 


0.0 


0.0 


0.0 


1.3 


1.3 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


Syst. 


5.7 


5.3 


4.9 


3.9 


3.6 


1.3 


2.9 


2.6 


3.0 


1.5 


1.7 


2.1 


Stat. 


10.0 


5.1 


10.3 


3.6 


12.3 


1.0 


2.7 


3.3 


6.2 


1.5 


1.3 


3.2 


Total 


11.5 


7.3 


11.4 


5.3 


12.8 


1.7 


4.0 


4.2 


6.9 


2.1 


2.2 


3.9 



Table 1: Summary of the measurements used to determine the world average M t . Integrated 
luminosity (/ C dt) is in fb _1 , and all other numbers are in GeV/c 2 . The error categories and 
their correlations are described in the text. The total systematic uncertainty and the total 
uncertainty are obtained by adding the relevant contributions in quadrature. 

determined from external calibration samples. Some parts of the associated uncertainty are 
correlated with the Run I JES uncertainty as noted below. 

The D0 Run Ha 1+j analysis is using the JES determined from the external calibration 
derived using 7+jets events as an additional Gaussian constraint to the in-situ calibration. 
Therefore the total resulting JES uncertainty has been split into the part coming solely from 
the in-situ calibration and the part coming from the external calibration. To do that, the 
measurement without external JES constraint has been combined iteratively with a pseudo- 
measurement using the BLUE method that would use only the external calibration so that the 
combination gives the actual total JES uncertainty. The splitting obtained in this way is used 
to assess the iJES and dJES uncertainties |29|. 
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A new analysis technique from CDF is included (trk) . This measurement uses both the mean 
decay-length from B-tagged jets and the mean lepton transverse momentum to determine the 
top-quark mass in 1+j candidate events. While the statistical sensitivity is not as good as 
the more traditional methods, this technique has the advantage that since it uses primarily 
tracking information, it is almost entirely independent of JES uncertainties. As the statitistics 
of this sample continue to grow, this method could offer a nice cross-check of the top-quark 
mass that's largely independent of the dominant JES systematic uncertainty which plagues the 
other measurements. The statistical correlation between an earlier version of the trk analysis 
and a traditional Run II CDF 1+j measurement was studied using Monte Carlo signal-plus- 
background pseudo-experiments which correctly account for the sample overlap and was found 
to be consistent with zero (to within < 1%) independent of the assumed top-quark mass. 

The two D0 Run II lepton+jets results [TTI [T8] are derived from Run Ha and Run lib 
datasets, respectively, and are labelled as such. The D0 Run II dilepton result is itself a 
combination of two results using different techniques but partially overlapping dilepton data 
sets [HIED]. 

Table [T] also lists the uncertainties of the results, sub-divided into the categories described 
in the next Section. The correlations between the inputs are described in Section HI 



3 Error Categories 

We employ the same error categories as used for the previous world average [28], plus one new 
category (lepPt). They include a detailed breakdown of the various sources of uncertainty and 
aim to lump together sources of systematic uncertainty that share the same or similar origin. For 
example, the "Signal" category discussed below includes the uncertainties from ISR, FSR, and 
PDF — all of which affect the modeling of the tt signal. Some systematic uncertainties have been 
broken down into multiple categories in order to accommodate specific types of correlations. 
For example, the jet energy scale (JES) uncertainty is sub-divided into several components in 
order to more accurately accommodate our best estimate of the relevant correlations. Each 
error category is discussed below. 

Statistical: The statistical uncertainty associated with the M t determination. 

iJES: That part of the JES uncertainty which originates from in-situ calibration procedures 
and is uncorrelated among the measurements. In the combination reported here it cor- 
responds to the statistical uncertainty associated with the JES determination using the 
W — > qq' invariant mass in the CDF Run II 1+j and all-h measurements and D0 Run Ha 
and Run lib 1+j measurements. Residual JES uncertainties, which arise from effects not 
considered in the in-situ calibration, are included in other categories. 
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aJES: That part of the JES uncertainty which originates from differences in detector e/h 
response between 6-jets and light-quark jets. It is specific to the D0 Run II measurements 
and is taken to be uncorrelated with the D0 Run I and CDF measurements. 

bJES: That part of the JES uncertainty which originates from uncertainties specific to the 
modeling of 6-jets and which is correlated across all measurements. For both CDF and 
D0 this includes uncertainties arising from variations in the semi-leptonic branching 
fraction, ^-fragmentation modeling, and differences in the color flow between 6-jets and 
light-quark jets. These were determined from Run II studies but back-propagated to the 
Run I measurements, whose rJES uncertainties (see below) were then corrected in order 
to keep the total JES uncertainty constant. 

cJES: That part of the JES uncertainty which originates from modeling uncertainties corre- 
lated across all measurements. Specifically it includes the modeling uncertainties associ- 
ated with light-quark fragmentation and out-of-cone corrections. 

dJES: That part of the JES uncertainty which originates from limitations in the calibration 
data samples used and which is correlated between measurements within the same data- 
taking period, such as Run I or Run II, but not between experiments. For CDF this 
corresponds to uncertainties associated with the 77-dependent JES corrections which are 
estimated using di-jet data events. For D0 this includes uncertainties in the calorimeter 
response to light-quark jets, and rj- and ^-dependent uncertainties constrained using 
Run II 7+jet data samples. 

rJES: The remaining part of the JES uncertainty which is correlated between all measurements 
of the same experiment independent of data-taking period, but is uncorrelated between 
experiments. For CDF, this is dominated by uncertainties in the calorimeter response 
to light-quark jets, and also includes small uncertainties associated with the multiple 
interaction and underlying event corrections. 

lepPt: The systematic uncertainty arising from uncertainties in the scale of lepton transverse 
momentum measurements. This is an important uncertainty for CDF's track-based mea- 
surement. It was not considered as a source of systematic uncertainty in the Run I 
measurements or in measurements at D0. 

Signal: The systematic uncertainty arising from uncertainties in the modeling of the it signal 
which is correlated across all measurements. This includes uncertainties from variations 
in the ISR, FSR, and PDF descriptions used to generate the it Monte Carlo samples 
that calibrate each method. It also includes small uncertainties associated with biases 
associated with the identification of 6-jets. 

Background: The systematic uncertainty arising from uncertainties in modeling the domi- 
nant background sources and correlated across all measurements in the same channel. 
These include uncertainties on the background composition and shape. In particular 
uncertainties associated with the modeling of the QCD multi-jet background (all-j and 
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1+j), uncertainties associated with the modeling of the Drell-Yan background (di-1), and 
uncertainties associated with variations of the fragmentation scale used to model W+jets 
background (all channels) are included. 

Fit: The systematic uncertainty arising from any source specific to a particular fit method, 
including the finite Monte Carlo statistics available to calibrate each method. 

Monte Carlo: The systematic uncertainty associated with variations of the physics model 
used to calibrate the fit methods and correlated across all measurements. For CDF it 
includes variations observed when substituting PYTHIA [3UJ ED E2] (Run I and Run II) 
or ISAJET [33] (Run I) for HERWIG [331 [35] when modeling the ti signal. Similar 
variations are included for the D0 Run I measurements. The D0 Run II measurements 
use ALPGEN [36J to model the it signal and the variations considered are included in 
the Signal category above. 

UN/MI: This is specific to D0 and includes the uncertainty arising from uranium noise in the 
D0 calorimeter and from the multiple interaction corrections to the JES. For D0 Run I 
these uncertainties were sizable, while for Run II, owing to the shorter integration time 
and in-situ JES determination, these uncertainties are negligible. 



These categories represent the current preliminary understanding of the various sources of 
uncertainty and their correlations. We expect these to evolve as we continue to probe each 
method's sensitivity to the various systematic sources with ever improving precision. Varia- 
tions in the assignment of uncertainties to the error categories, in the back-propagation of the 
bJES uncertainties to Run I measurements, in the approximations made to symmetrize the 
uncertainties used in the combination, and in the assumed magnitude of the correlations all 
negligibly effect (^C 0.1GeV/c 2 ) the combined M t and total uncertainty. 



4 Correlations 

The following correlations are used when making the combination: 

• The uncertainties in the Statistical, Fit, and iJES categories are taken to be uncorrected 
among the measurements. 

• The uncertainties in the aJES and dJES categories are taken to be 100% correlated among 
all Run I and all Run II measurements on the same experiment, but uncorrelated between 
Run I and Run II and uncorrelated between the experiments. 

• The uncertainties in the rJES and UN/MI categories are taken to be 100% correlated 
among all measurements on the same experiment. 
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Run I published 


Run II preliminary 


CDF 
1+j di-1 all-j 


D0 
1+j di-1 


CDF 

1+j di-1 all-j trk 


D0 

1+j/a 1+j/b di-1 
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CDF-I di-1 
CDF-I aU-j 


1.00 

0.29 1.00 

0.32 0.19 1.00 
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0.26 0.15 0.14 
0.11 0.08 0.07 


1.00 

0.16 1.00 






CDF-II 1+j 
CDF-II di-1 
CDF-II aU-j 
CDF-II trk 


0.32 0.18 0.20 
0.45 0.28 0.33 
0.17 0.11 0.16 
0.16 0.08 0.07 


0.19 0.07 
0.22 0.11 
0.10 0.05 
0.13 0.05 


1.00 

0.34 1.00 

0.15 0.19 1.00 

0.17 0.12 0.05 1.00 




D0-II 1+j/a 
D0-II 1+j/b 
D0-II di-1 


0.11 0.05 0.03 
0.12 0.06 0.04 
0.05 0.04 0.02 


0.08 0.03 
0.09 0.03 
0.04 0.04 


0.09 0.04 0.03 0.09 
0.10 0.05 0.03 0.09 
0.04 0.05 0.03 0.03 


1.00 

0.24 1.00 

0.31 0.16 1.00 



Table 2: The resulting matrix of total correlation coefficients used to determined the world 
average top quark mass. 

• The uncertainties in the Background category are taken to be 100% correlated among all 
measurements in the same channel. 



• The uncertainties in the bJES, cJES, Signal, and Generator categories are taken to be 
100% correlated among all measurements. 



Using the inputs from Table [T] and the correlations specified here, the resulting matrix of total 
correlation co-efficients is given in Table [2j 

The measurements are combined using a program implementing a numerical y 2 minimiza- 
tion as well as the analytic BLUE method [21], [22]. The two methods used are mathematically 
equivalent, and are also equivalent to the method used in an older combination [37], and give 
identical results for the combination. In addition, the BLUE method yields the decomposition 
of the error on the average in terms of the error categories specified for the input measure- 
ments [22]. 
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Run I published 


Run II preliminary 


CDF 
1+j di-1 all-j 


D0 
1+j di-1 


CDF 

1+j di-1 all-j trk 


D0 

1+j/a 1+j/b di-1 


Pull 

Weight [%} 


+0.5 -0.4 +1.2 
-3.4 -0.6 -0.6 


+1.5 -0.3 
+1.2 +0.2 


-0.1 -0.3 +1.1 +0.4 
+46.1 +3.7 +5.2 +0.02 


-0.5 +0.3 +0.5 
+23.7 +21.5 +3.0 



Table 3: The pull and weight for each of the inputs used to determine the world average mass 
of the top quark. See Reference [2T] for a discussion of negative weights. 

5 Results 

The combined value for the top-quark mass is: 

M t = 172.4 ± 1.2 GeV/c 2 , (1) 

with a x 2 °f 6.9 for 11 degrees of freedom, which corresponds to a probability of 81%, indicating 
good agreement among all the input measurements. The total uncertainty can be sub-divided 
into the contributions from the various error categories as: Statistical (±0.7), total JES (±0.8), 
Lepton scale (±0.1), Signal (±0.3), Background (±0.3), Fit (±0.1), Monte Carlo (±0.3), and 
UN/MI (±0.02), for a total Systematic (±1.0), where all numbers are in units of GeV/c 2 . The 
pull and weight for each of the inputs are listed in Table [3j The input measurements and the 
resulting world average mass of the top quark are summarized in Figure [TJ 

The weights of many of the Run I measurements are negative. In general, this situation 
can occur if the correlation between two measurements is larger than the ratio of their total 
uncertainties. This is indeed the case here. In these instances the less precise measurement 
will usually acquire a negative weight. While a weight of zero means that a particular input 
is effectively ignored in the combination, a negative weight means that it affects the resulting 
central value and helps reduce the total uncertainty. See reference [21] for further discussion of 
negative weights. 

Although the x 2 from the combination of all measurements indicates that there is good 
agreement among them, and no input has an anomalously large pull, it is still interesting to 
also fit for the top-quark mass in the all-j, 1+j, and di-1 channels separately. We use the same 
methodology, inputs, error categories, and correlations as described above, but fit for the three 
physical observables, Mf~ l ~\ Ml + \ and M t dl_1 . The results of this combination are shown in 
Table H] and have x 2 °f 4.9 for 9 degrees of freedom, which corresponds to a probability of 84%. 
These results differ from a naive combination, where only the measurements in a given channel 
contribute to the M t determination in that channel, since the combination here fully accounts 
for all correlations, including those which cross-correlate the different channels. Using the 
results of Table @] we calculate the chi-squared consistency between any two channels, including 
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Mass of the Top Quark (*Preliminary) 



CDF-I di-l 



DO-I di-l 
CDF-II di-l 
*D0-ll di-l 
CDF-I l+j 
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DO-llb l+j 
CDF-I all-j 
CDF-II all-j 
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Tevatron July'08 




- 



167.4 ±10.3 ±4.9 
168.4 ±12.3 ±3.6 
171 .2 ±2.7 ±2.9 
174.4 ±3.2 ±2.1 



176.1± 5.1 ± 5.3 
180.1±3.6±3.9 
172.2± 1.0± 1.3 
1 71 .5 ± 1.5 ± 1.5 
1 73.0 ± 1.3 ± 1.7 



186.0 ±10.0 ±5.7 
176.9 ±3.3 ±2.6 

175.3 ±6.2 ±3.0 

1 72.4 ± 0.7 ± 1.0 

(stat.) ± (syst.) 
X 2 /dof = 6.9/1 1 .0 (81%) 
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Figure 1: A summary of the input measurements and resulting world average mass of the top 
quark. 
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Parameter 


Value (GeV/c 2 ) 


Correlations 


M all-j 

M t 1+j 
M? 1 - 1 


177.5 ±4.0 
172.2 ± 1.2 
171.5 ±2.6 


1.00 

0.14 1.00 

0.17 0.32 1.00 



Table 4: Summary of the combination of the 12 measurements by CDF and D0 in terms of 
three physical quantities, the mass of the top quark in the all-jets, lepton±jets, and di-lepton 
channels. 

all correlations, as x 2 (dil — Ij) = 0.08, x 2 (^J ~ a>Uj) = 1-7, and x 2 ( a ^j ~ dil) = 1.9. These 
correspond to chi-squared probabilities of 78%, 19%, and 17%, respectively, and indicate that 
the determinations of M t from the three channels are consistent with one another. 



6 Summary 

A preliminary combination of measurements of the mass of the top quark from the Tevatron 
experiments CDF and D0 is presented. The combination includes five published Run I mea- 
surements and seven preliminary Run II measurements. Taking into account the statistical 
and systematic uncertainties and their correlations, the preliminary world-average result is: 
M t = 172.4 ± 1.2 GeV/c 2 , where the total uncertainty is obtained assuming Gaussian sys- 
tematic uncertainties and adding them plus the statistical uncertainty in quadrature. While 
the central value is somewhat higher than our 2007 average, the averages are compatible as 
appreciably more luminosity and refined analysis techniques are now used. 

The mass of the top quark is now known with a relative precision of 0.7%, limited by 
the systematic uncertainties, which are dominated by the jet energy scale uncertainty. This 
systematic is expected to improve as larger data sets are collected since new analysis techniques 
constrain the jet energy scale using in-situ W — * qq' decays. It can be reasonably expected 
that with the full Run II data set the top-quark mass will be known to better than 0.7%. 
To reach this level of precision further work is required to determine more accurately the 
various correlations present, and to understand more precisely the 6-jet modeling, Signal, and 
Background uncertainties which may limit the sensitivity at larger data sets. Limitations of 
the Monte Carlo generators used to calibrate each fit method also become more important as 
the precision reaches the ~ 1 GeV/c 2 level; these warrant further study in the near future. 
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